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Abstract 

This paper presents innovative algorithms to efficiently compute erosions and dilations of run-length encoded (RLE) 
binary images with arbitrary shaped structuring elements. An RLE image is given by a set of runs, where a run 
is a horizontal concatenation of foreground pixels. The proposed algorithms extract the skeleton of the structuring 
element and build distance tables of the input image, which are storing the distance to the next background pixel on 
the left and right hand sides. This information is then used to speed up the calculations of the erosion and dilation 
operator by enabling the use of techniques which allow to skip the analysis of certain pixels whenever a hit or miss 
occurs. Additionally the input image gets trimmed during the preprocessing steps on the base of two primitive 
criteria. Experimental results show the advantages over other algorithms. The source code of our algorithms is 
available in C++. 

Keywords: binary image, dilation, erosion, filtering algorithms, image analysis, image denoising, computer vision, 
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1. Introduction 

Mathematical morphology is a general method for the analysis of spatial structures which aims at 
analysing the shape and form of objects [1]. In a variety of industrial computer vision applications, ranging 
from barcode scanning to the placement of chips in semiconductor industry, mathematical morphology is 
being used to process images and filter noise. We are interested in analysing binary images since these 
can be represented as sets. Therefore it enables us to use set-theoretical tools to process these images. 
Morphological operators which are used for noise filtering are constructed using two basic operators, namely 
erosion and dilation. In this paper we propose fast erosion and dilation algorithms on two-dimensional 
run-length encoded (RLE) binary images. A run R = {lx, rx, y) is a concatenation of pixels in horizontal 
direction, where y indicates the y—coordinate of these pixels, lx the a;—coordinate of the leftmost pixel 
and rx the a;-coordinate of the rightmost pixel, formally R = {lx, rx, y) = {{lx, y), {lx l,y),..., {rx, y)}. 
Then a two-dimensional binary image Z can also be described by its run-length representation 


Previously, several very different ideas for the fast computation of erosion and dilation have been presented, 
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such as using the decomposition of a rectangular-shaped structuring element [2], basing the algorithm on 
set-theoretical investigations [3, 4], solely operating on the contours of the image [5], or applying methods 
that are similar to string-matching techniques [6]. Inspired by some of the ideas introduced in [3, 4, 6] we 
propose two algorithms that extend these ideas by a couple of new theorems. 

An outline of our paper is as follows: In Section 3 we provide the mathematical investigation of the erosion 
operator on RLE images. We come up with various ideas which can be used to formulate a fast algorithm. 
In Section 4 we take use of the duality between erosion and dilation to port our ideas to not only work 
with erosion but also with dilation. In Section 5 we propose fast algorithms based on the theorems given 
in the previous sections and in Section 6 we analyse the complexity of these algorithms. Finally we present 
runtime experiments and comparisons with other algorithms in Section 7 and our conclusion can be found 
in Section 8. 

2. Preliminaries 

As described in the Introduction, morphological operators are constructed using the two basic operators: 
erosion and dilation. A binary image is given by a set of pixels. Assuming X C and the structuring 
element (SE) B C to be binary images, these operators are defined as: 


£b{X) = 

[peZ^lBpCxy 

(1) 

Sb{X) = 

{pGZ‘'|(R‘)pnA^0}, 

(2) 


where £b{X) denotes the erosion of X by B, Sb{X) the dilation of X by B, B* the reflection about the 
origin and Bp the translation of i? by a vector p G Instead of eB(A) and Sb{X) we also write XqB and 
X (B B whenever it is convenient. A two-dimensional binary image Z can also be described by its run-length 
representation Z = where N denotes the number of runs and i?„ the nth run of Z. A compact 

representation of an RLE image is given when concatenated pixels are defined as a single run (not divided 
into several runs) and when runs are not overlapping. The proposed algorithms require compact RLE binary 
images as input and also return the eroded or dilated image in compact RLE representation. 

3. Erosion Using RLE 

In this chapter, we develop a theory to speed up the calculation of the erosion operator. This is done 

in five steps. In Section 3.1 we formulate a theorem which allows us to translate the structuring element 

before calculating the erosion such that the translated structuring element contains the origin and describe, 

how this improves the efficiency of our algorithm. In Section 3.2 we formulate a primitive criterion which 
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allows us to remove all runs of the input image that are shorter than the longest run within the structuring 
element. This is achieved in linear runtime complexity before calculating the erosion. Section 3.3 on one 
hand summarizes the results of [6] which are needed to prove the lemmata and theorems stated in Sections 
3.4 and 3.5, on the other hand it also describes an efficient way to build the erosion transform tables and 
skeletons of a given image. In Section 3.4 we extend the Jump-Miss Theorem by formulating the Jump-Hit 
Theorem and in Section 3.5 we give another - easy to apply - criterion, which makes the investigation of a 
big number of pixels redundant (depending on the length of the shortest run of the structuring element). 

3.1. Translating the SE such that it contains the origin 

When the origin o = (0,..., 0) G is contained in H (o G B) expression (1) obviously reduces to 

eB(X) = {pGX|HpCx}. (3) 

The following theorem allows us to translate H by a vector q such that o G Bg before calculating the erosion. 
Theorem 1. Let X, B C and q . Then we have that: 

£b{X) = [£B^{X)]q. 


Proof. Obviously 


[£Bg{X)]q — {p I Bp+q C X}q 

= {p + q\ Bp+q C X} 

= {p I Bp ^ X} by substitution p \= p -\- q 

= £b{X). 


Since our image X is given in RLE representation, we are able to look up d\\ p G X efficiently. 

3.2. Removing short runs before eroding 

Let X = Rn and B = lJm=i be two run-length encoded binary images in compact repre¬ 

sentation. The article [3] provided an elegant method to calculate the erosion directly using the RLE 
representations of B and X: 

M N 

£BiX)= f] \JenB{R^). (4) 

m—1 n—1 

While the erosion of runs can be computed very efficiently, the intersections and unions of runs take more 
effort. A runtime comparison with an algorithm which is based on this method can be found in Section 7. 
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We are able to use this result to prove the following theorem. It allows us to remove all runs of X that are 
shorter than the shortest run within B before calculating the erosion. 

Theorem 2. Let X = ^ — Um=i binary images in compact RLE representation 

and 


Lmin '■= min ll-Rmll as Well as 
Xl . := 11 R^. 

^min El 

l<n<N 

Then we have that: 

£b{X) = £b{Xl^.^). 

Proof. For all runs with \Rn\ < Tmin and all m € {1,..., M} the following equality obviously holds: 


^R^{Rn) — 0- 


( 5 ) 


In the following, let if” := eii^{R^)- By applying (4) we get: 


M N 


^ b { x ) = f| IJ if" 


m—1 n—1 
M / r 

= n UK 

m—1 \ l<n<A^ 

M 

= n U e : 

m—1 l<n<A^ 


u 


U K 


l<n<Af 


using (5) 


= £B{XL^i„)- 


3.3. Jump-Miss Theorem 


Let X.AQ'E’^ be two binary images with |ff| < oo, and |A| > 1. 


Then the n-fold erosion is given by 


X 0„ A 


j X for n = 0, 

(X Qn-i A) Q A for n > 0. 


4 






G. Ehrensperger ■ A. Ostermann ■ F. Schwitzer / 00 (2015) 1-17 


5 


This allows us to introduce the erosion transform by 

{ max {n \ p e X 0„_i A} if p £ X, 

n>0 

0 otherwise. 

The skeleton Sg of B with respect to A (where |A| > 1 and o G A) is the set of all points p G B which 
satisfy the following inequality (see [6, Proposition 2.7]): 

max {f^{p + e)} < f^{p). (6) 

eG A* 

According to the following theorem when checking whether a structuring element B is contained in A at a 
position G it is enough to compare the values /b^(s + h) and fx{s + h) at all points of the skeleton 
sG 5^. 

Theorem 3. [6, Theorem 3.3] Let A,B,X C |A|, \B\, |A| < oo with o G A, |A| > 1, and h G Then 
the following are equivalent: 

1. G £b{X') 

2. Bh<LX 

3. \/s G Sg : fs^ih + s) < fxih + s). 

In practice, for most structuring elements B we have that jS”^] <C |i?|. So once the erosion transform 
and skeleton are computed, the runtime of the actual erosion reduces to (!1(|A||S'^|). Furthermore it was 
shown that it is also possible to skip the analysis of certain pixels whenever fs^ih + s) — fx{h + s) > 0. 
This fact is described by the Jump-Miss Theorem which will be discussed next. 

Theorem 4. (Jump-Miss Theorem) [6, Theorem 3.5] Let A,B,X C with o G A and |A|, |i?|, |A| < 
oo, |A| > 1 as well as s G and h G TA . Then 

fx{h + s) < /b,^(^ + s) 


implies that there are at least 

k ■■= fsi^ih + s) - f^{h + s) 


points g G TA such that Bg ^ X. These k points are given by gi := h i ■ e for all i G {0,..., fc — 1} and 
e G A*. 


Now the question arises whether the erosion transform f^^ needs to be recomputed for every given h. 

The following lemma states that it is sufficient to only determine fg. 
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Lemma 5. [6, Proposition 3.2] Let A, B G with \A\ > 1, o G A and |A|, \B\ < oo. Then for all h G 
and s G Sg: 

/ b ( s ) = fs^h + s). 

When thinking of a convenient set A to work with on run-length encoded images, one immediately comes 
up with A = {(—1,0), (0,0)} or a similarly shaped form. This set ensures fast computation of /xi/s ^-nd 
allows one to determine Sg very efficiently. With this set the erosion transform evaluates to 0 outside of 
X, to 1 for the left-most pixel of a given run R]] C X and increments its value for every following pixel (in 
horizontal direction) within the run, such that evaluates to | for the rightmost pixel. The following 
theorem proves this statement. 

Theorem 6. Let A,XG1 |X| < oo, h G X, and A = {(—1, 0), (0,0)}. Then 

fx{h) = fxih + (—1,0)) + 1. 


Proof. We distinguish between two cases: 

1. Case 1 {h + (—1,0) ^ X): 

Then by definition we have that /^(/i -|- (—1, 0)) =0, h G X = X Qq A and obviously 

xeiT = {p|{(-i,o),(o,o)}pCx} 

= {p I {(- 1 , 0 ) (0,0) + p} x) ^ h 

because by assumption ((—1, 0) + h) ^ X. Thus fx{h) = 1 = fx{h + (—1, 0)) -I- 1. 

2. Case 2 (/i-H (-1,0) G X): 

Without loss of generality let fx{h + {—1,0)) = h. First we prove the following implication by 
induction: 


{hGX and h+ (-1,0) G X Qn-i A) h G X Qa A. (7) 

(a) Base case (n = 1): When h,h+ (—1,0) G X, obviously also hGXQA = {p\ {(—1, 0), (0,0)}p C 

x} = {p\{p + {-i,0),p}cx). 

(b) Induction step (n -|- 1): Let h G X and h + (—1,0) G X Qn A. Since o G A, erosion is an anti- 
extensive operator [1], thus also /i-|- (—1,0) G X Qfi-i A and by induction hypothesis h G XQ^ A. 
We conclude X Qn+x A = ]p \ {{p + {—1, 0),_p)} (G X Qn A] ^ h. This proves (7). 

From fx{h + (—1,0)) = n we conclude that h + (—1,0) G X A and h + (—1, 0) ^ X 0^^ A. Thus 

0n+l A= ]p\ {p + (—1, 0),p} C X Qn ^ h 
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and using (7), we get h £ X Qn A. Putting all together, this implies fx{h) = h + 1. 

Note that (6) and Theorem 6 imply that the skeleton can be retrieved by taking the rightmost pixel of 
each run. 

3.4- Jump-Hit Theorem 

In the previous section we described a criterion which allows us to skip the analysis of certain pixels 
whenever the algorithm finds a pixel h such that h ^ £b{X). We also call this a miss. Now we formulate a 

criterion which applies for hits. A hit denotes the occurrence of an h such that h € eg (A"). 

Lemma 7. Let X C he a binary image, |X| < oo, p G X and A = {(—1, 0), (0, 0)}. Then the condition 
fx (p) — ^ implies 

Vi G {0,..., n - 1} : fl4{p + (i, 0)) = f^{p) + i. 

Proof. Obviously this statement holds for i = 0. Let i G {1,.. ., n — Ij. For the given set A the equality 
fx\p) — implies that for all z G {1,..., n — 1} also (p + (i,0)) G X. By applying Theorem 6 multiple 
times we get f^{p + (i, 0)) = /^(p) + i. 

Note that Lemma 7 does not hold for arbitrary sets A. 

Theorem 8. (Jump-Hit Theorem) Let B,X C 1?, A = {(—1,0), (0, 0)} and {hx,hy) = h G 'I?. If for 

all sgS^ we have that fg^^ih + s) < fxi^ + s), then the complete eroded run is given by 


{hx: l^x T n 1, hy'^ C sB(,Xf 


where n 


ram{f^\h + s)}. 


Proof. Using Theorem 3 we get h G eB(X) and Lemma 7 implies for all s G and all i G {0,..., n — 1} 
that 


fx{h + S + (i, 0)) = fx{h + s) + Z with Lemma 7 

>fxih + s) i>0 

— /sfe T s), by assumption 

from which we conclude - again with Theorem 3 - that {h + (z, 0)) G e_B(X) for all z G {0,..., rz — 1} and 
therefore {hx, hx -G n — 1, hy) C eB(X). 

We want to emphasize that, in case of a miss, the length of a jump is limited by the width of the 

structuring element. When having a hit, however, the length of the jump is always maximal and determined 

by just checking the values of fx*{h + s) at its starting pixel h along the pixels s of the skeleton Sg. 

7 


G. Ehrensperger • A. Ostermann • F. Schwitzer / 00 (2015) 1-17 


3.5. Skipping the analysis of additional pixels 

The goal of this section is to find a binary image Xcut C such that |Xcut| ^ and 


In the upcoming lemma we investigate a very special class of structuring elements. These elements 
include a run whose rightmost element is the origin. In a next step we are going to extend the lemma such 
that this condition is no longer needed. 

Lemma 9. Let X = lj^=i Rn R = Um=i R-m run-length encoded binary images in compact 

representation, ifi & {1,..., M} such that i?? = {Ixm, 0, 0), and X := , where 

~X \{^Xn + {\R^\-l),rXn,yu) iI\R^\<\Rn\, 

■ \ 

I 0 else 

for all n = 1,.. ., N. Then we have that: 

SB{X) = [pex\BpCXL„,,„}, 


where defined as in Theorem 2. 

Proof. Let be an arbitrary run of X, A := {(—I, 0), (0, 0)}, and h G R^. Then, since R^ = {Ixm, 0, 0), 
using Theorems 1 and 6, the following equality holds: 


flih+i0,0)) = f^,^{h) = /i((0,0)) = \Rg\. 


We distinguish between two cases: 

1. Case I: Let \Rg\ > \Rg\- Then, by using the previous equality and Theorem 6, we get: 

fl{h + {0,0)) = \Rg\>\Rg\>f^{h). 

Theorem 3 implies h ^ esiX) for all /i G Rg . 

2. Case 2: Let \Rg\ < \Rg\- Using Theorem 6 we conclude that for the first \Rg\ — 1 pixels h of run Rg 
we have that fx{R) < \Rg,\- This leads to: 

fl{h + {0,0)) = /^((0,0)) = \Rg\ > f^{h) 

for/i G {lxn,lxn + \Rg\ — ‘2,yn)- Again Theorem 3 implies/i ^ eB(X) for all/i G {lxn,lxn + \Rg\ — ‘^,yn)- 
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Since o £ B hy assumption, (3) and Theorem 2 state that esiX) = {p G \ Bp C Additionally 

we just proved that for all n S N} and h £ {Ixn, Ixn + \R^\ — 2, ?/„) where (a, b, c) := 0 for a > 5 also 

h ^ £b{X). Thus the statement holds. 

Obviously it makes sense to translate B by a vector q such that the rightmost pixel of the longest run 
within B is placed at the origin. This translation is described in the next theorem. 

Theorem 10. Let X = and B = R^ be two run-length encoded binary images in compact 

representation and Aleut := where max {|.Rml} 

~X \ + (Ajuax ~ if < \Rj^ I, 

Rji ^ 

I 0 else 

for all n = 1,. ..,N. Let i?? C B be any run such that \R%\ = Tmax- Then q := (ra:?,y?) denotes the 
rightmost pixel of the longest run R^ within B and we have that: 


SBiX) = [eb.M)]-, = {p G ^cut I Bp_, C 


where Xb^^^ is defined as in Theorem 2. 


Proof. Theorem 1 allows us to translate the structuring element B hy q = {rx^,yi^). So o G B-q and 
{R^)-q = {Ix^ — rx^,0,0). Hence the claim immediately follows by using Lemma 9. 


Altough we can reduce our investigation to the points of Aleut, we still need to evaluate and 

fxj^ . Figure 1 demonstrates the efficiency of this theorem. 


4. Dilation Using RLE 

By taking use of the duality between erosion and dilation [1] 


Sb{X) = 


c 


SB^X^) , 


we are able to propose a fast dilation algorithm based on the previously developed erosion algorithm in a 
straightforward way. But we have to make some additional thoughts, because obviously |A1'^| = oo when 
\X\ < oo. 

Consider the rectangle REC*’'^’*’^ with l^r^t,b £ 'L given by: 


■.= {{x,v)£'I?\l<x<r, t<y<b}. 
9 
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Then we denote the smallest rectangle which includes X by REC™“. 

Obviously (by considering the definition of dilation (2)) there exists a sufficiently big rectangle REC^ 
such that SsiX) does not contain any points outside this rectangle: 


6BiX)\RE& 




\ REC^ = 0. 


Therefore we get: 




n REC^ = eBt{X^) 


Because of | REC^ | < oo we are also allowed to restrict by a rectangle REC”^, such that 


eB*(^"nREC'^) nREC^= . 


( 8 ) 


Still we have to find snitable sizes for the rectangles REO^ and REC^. Our goal is to make them as small 
as possible while being able to compute their sizes efficiently. 

Because of Theorem 1 we can assume that o G B. If this is not the case, we are allowed to translate R by a 
suitable vector q such that o G Bq. Since 


Sb{X) = {pg^^ I {B%nx^n 


one can easily see, that a suitable rectangle REC^ is received by adding a border to REG™" which has the 
same width as the structuring element on the left and right hand sides and the same height as the structuring 
element on the top and bottom. Now we only need to find a rectangle REG® such that (8) holds. Therefore 
we just double the border of REG^ to obtain REG”^. 

With these thoughts we are able to reformulate the dilation algorithm using the proposed erosion algorithm: 


1. We receive run-length encoded images X and B as input. 

2. We calculate R* and the compact run-length representation of X‘^ G REG^. 

3. These images can be used to retrieve G REG^) by executing GenErosionTransX2cut(X'^ G 

REG^R‘). 


4. 


In the last step we take the complement of the erosion and restrict the image to the size of the rectangle 
REG^: 


Sb{X) 


eBt(X^GREG^) 


C 


G REG"^. 


Please note that 2) and 4) can also be done on the fly as described in Section 5.2. 
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Algorithm 1 GenerateSkeletonB 


1 

function GenerateSkeletonB(B) 

t> (01) 

2 



3 



4 

for all C B do 

> (02) 

5 

•S's ■*- U {(ra®,y^)} 

t> (03) 

6 

^ I-Rmi 

> (04) 

7 


t> (05) 

8 

end for 


9 

return (List(S^), /^, 

t> (06) 

10 

end function 



comments: 



(01): returns List(5^), fg and LA^ 

(02): is defined as = (Zx® , rr® , 

(03): adds the rightmost element to Sg 

(04): erosion transform of this element equals the length of the 



run 



(05): keeps track of the shortest run within B 
(06): List (5^) means that is stored as a list 



5. Proposed Algorithms 

In this section we are going to describe fast erosion and dilation algorithms. Reference [7] provides 
several figures that demonstrate the principles of these algorithms. 


5.1. Proposed Erosion Algorithm 

Observation (3) and Theorem 1 allow us to focus our investigations at those pixels of A, which can 
be looked up very efficiently for RLE images. Due to Theorem 2, we are allowed to drop all runs of X 
during preprocessing which are shorter than the shortest run of B. Theorem 10 states that it is enough to 
investigate the pixels of Acut, which is obtained by removing the first ~ 1 pixels of every run of A, 

where L^ax denotes the length of the longest run of B. Of course, runs that are shorter than vanish 

completely. In case of a miss, the Jump-Miss Theorem enables us to skip the analysis of certain pixels and 
whenever a hit occurs, the full eroded run can immediately be added by applying the Jump-Hit Theorem. 

A pseudocode that implements all of these methods is given by the Algorithms 1, 2, and 3. 


5.2. Proposed Dilation Algorithm 

By adapting the described erosion algorithm, all required operations (calculating the complement, inter¬ 
secting with rectangles) can be done on the fly. Recall that we obtained the following formula in Section 4: 


5b{X)= £Bt(A^nREC'') 


n REC"^. 
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Algorithm 2 GenErosionTransX2cut 

1 

function GENEROSIONTRANSX2CUT(Jf, imax) 

> (01) 

2 

for all 2 : G Mx do 

t> (02) 

3 



4 

/£ . 


5 

end for 


6 

-^cut ^ 0 


7 

for all C X do 

t> (03) 

8 

if \R^\ > then 

t> (04) 

9 

Xcut ^ Xcut U (Ix^ + - l,rx^,y^) 

> (05) 

10 

i 1 


11 

for i <— to rx^ do 

t> (06) 

12 

^XL^.X^Vn)) ■*- 1 


13 

fxL . ((hiln )) <- l-Rn 1 - i + 1 


14 

j <-7+1 


15 

end for 


16 

end if 


17 

end for 


18 

return , /jt* ,Xcut) 


19 

end function 



comments: 



(01): returns , /^‘ and JsTcut 

■^min -^min 

(02): initializes and ; Mx denotes a 

■^min -^min 

k X 1 array where k := width(X) and 1 := height(X) 



(03): visits every run in X; is defined as R(^ = 



{lx^,rx^,y^} 



(04): only considers runs with length > 

(05): where (a, b, c) := 0 when a > b 

(06): visits every point within the given run and calculates the 



according erosion transform values 



By replacing line 5 in GenerateSkeletonB with ^ U {(—“l/m)} ^ble to generate the 

skeleton of B* on the fly. In GenErosionTransX2cut we immediately generate the erosion transform of 
by adding the sequences between two runs C X instead of the runs themselves. Additionally 
we need to add the border with width and height of B. Last but not least, we can modify GetErosion2cut 
by adding the misses instead of the hits to the dilated image. 


6. Runtime Analysis 

6.1. Erosion 

In Section 3.3 we stated that the runtime complexity of the actual erosion algorithm is (!l(|X||S'g|). Next 
we are going to investigate the preprocessing steps. To build the skeleton Sg we need to visit every run 
and store the length of the run as well as the rightmost pixel. Then the complexity is given by 0{M), 
where M denotes the number of runs of B and obviously this leads to |iS'^| = M. To obtain and 
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Algorithm 3 GetErosion2cut 


1 

function GetErosion2cut(JC, B) 

> (01) 

(01): returns £b{-^) 

2 

I'max 0 


(02): determines the length and the position 

3 

for all C B do 

> (02) 

of the rightmost pixel of the longest 

4 

if f-max < |f?ml then 


run within B 

5 


> (03) 

(03): sets q to the rightmost pixel of the 

6 

7 

q {rx^,y^) 

end if 

longest run 

8 



(04): B—q contains the origin, which is the 

9 

B-q ■<— TranslateImage(S, —q) 

> (04) 

rightmost pixel of the longest run of B 

10 

List(54 ),/b , <—GenerateSkeletonB(B_q) 


(05): visits every run of Acut; Rn defined 

11 

fxr ’fxl , Xcut <-GenErosionTransX2cut(X, LA 

as = {lx(( ,rx^ ,y^) 

12 

min min 


(06): stores the x-coordinate of point h\ we 

13 

for all C Xcut do 

> (05) 

are checking if B^—q Q A 

14 

X ^ Ix^ 

> (06) 

(07): breaks the loop once our x lies outside 

15 

while X < rx^ do 

t> (07) 

the given sequence 

16 

miss <— false 

t> (08) 

(08): variable needed for breaking the loop 

17 

s •<— head(List(Sg )) 

> (09) 

in case of a miss according to the 

18 

while ((not(miss)) AND (exists(next(List(5^ ))))) do 

t> (10) 

Jump-Miss Theorem 

19 

s -i— next(List(5^ )) 

t> (11) 

(09): s points at the head of List(5R ) 

20 

Dlff ^ /s (s) - /x. + 

5 min 


(10): loop will be continued as long as there 

21 

while ((Dili > U) AND [x < rx^)) do 

t> (12) 

is no miss and there are more s G 54 

22 

23 

X X Diff 

Diff ^/s (s) -/x^ (s + {x,y(()) 

> (13) 

B-q 

which have not been visited yet 

24 

miss ^ true 

t> (14) 

(11): gets the next element within 

25 

end while 

List(Sg ) 

26 

end while 


(12): 

27 

if (not(miss)) then 

t> (15) 

(13): jump by Diff number of pixels 

28 

minDist ^ oo 


(14): breaks the outer loop since B^—q ^ A 

29 

for all s G S'd do 
^-g 

> (14) 

(15): Bh—q C A, so the Jump-Hit Theorem 

30 

minDist ^ min(minDist,(s + {x,y^))) 

> (16) 

can be applied 

31 

end for 


(16): determines the minimum distance to 

32 

<- £s__j(X) U (x, X + minDist - l,y(() 

> (17) 

the right to be able to apply the 

33 

34 

35 

X X minDist + 1 

end if 

end while 

> (18) 

Jump-Hit Theorem 
(17): Jump-Hit Theorem 

36 

end for 


(18): next point (h — q) to look at will be at 

37 

£:r(A)T ranslatelmagef^R o) 

t> (19) 

a;-coordinate (x -|- minDist H- 1) 

38 

return (£b(X)) 


(19): we retrieve £:g(A) by applying 

39 

end function 


Theorem 1 


Ixl need to visit every point of X, thus its complexity is given by 0(|A|). Further Acut is built while 

determining and and does not affect the given complexity. Putting all together, we get a 

■^min -^min 

complexity of 


0{\X\M + M + \X\) = 0{\X\M). 


(9) 


13 







G. Ehrensperger • A. Ostermann • F. Schwitzer / 00 (2015) 1-17 


14 


Table 1. Abbreviations of the tested algorithms 


algorithm 

source 

RLE 

erosion as in [4] 

Machado 

erosion as in [6] 

OpenCV 

erosion as in [8] 

eEJMH 

erosion algorithm as in Sec. 5.1 

dRLE 

dilation as in [4] 

dOpenCV 

dilation as in [8] 

dEJMH 

dilation algorithm as in Sec. 5.2 


6.2. Dilation 

Since this algorithm is taking use of the duality between erosion and dilation, and building the comple¬ 
ments is done on the fly, we get the following complexity (compare with (9)): 

OdV^nREC'^l • M). 


7. Experimental Results 

In this section we are going to compare the runtimes of the proposed algorithms. Implementations of 
these algorithms are available in C++ and can be found at h.ttps://numerical-analysis.uibk.ac.at/g. 
ehrensperger. We are comparing them with the implementations of the free library OpenCV (developed by 
Intel) v2.4.2 and with the algorithms proposed in [4] and [6] . The latter two got implemented by Machado and 
the source code is available at http://score.ime.usp.br/~dandy/inestrado.php. This implementations 
are also used in our tests. Additionally he implemented various variants of his own algorithm. In the 
following plots we took the pointwise minimum of the runtimes of his algorithms. The conversion of X and 
B into the used input format of the various algorithms is not part of the given runtimes.^ 

All tests were carried out on a workstation with Intel Core i7-3770K (for the following tests only one core 
was used), 32 GB DDR3-1333 RAM and OS Ubuntu 13.04 64 Bit. CPU-stepping, overclocking settings, 
and various energy saving options had been disabled. The test environment was compiled with GCC (the 
GNU Gompiler Collection) v4.7.3 (official project URL: http://gcc.gnu.org/) and compiler flag 02. The 
following results are the arithmetic mean values of three iterations. We used the image in Figure 1 to 
compare the runtimes of the algorithms listed in Table 1. 

The results are given in Figure 2. Note that the plots’ y—axes are scaled logarithmically. As can be 
seen, the algorithms proposed in this paper tremendously improve the runtime of both, the erosion and 

^In [6] the conversion of the raster graphics X and B into the used input format of the algorithms was part of the given 
runtimes. 
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Figure 1. The size of this image is 2500 X 2500 pixels. The gray square (green in the online version) in the upper-left corner 
represents the structuring element, the gray parts (orange in the online version) of the image demonstrate the efficiency of 
Theorem 10. This theorem states that the investigation of the gray pixels can be skipped while calculating the erosion. The 
source of the original image is http://www.imageprocessingplace.com/. 

dilation operator, over the compared algorithms on the test image given in Figure 1. We observe that the 
execution times of RLE, dRLE, OpenCV, and dOpenCV grow with the size of the structuring element. 
In contrary eEJMH and Machado benefit from the size of the structuring element. This can be explained, 
since the bigger the structuring element the larger the jumps allowed by the Jump-Miss Theorem. The 
clear advantage of eEJMH over Machado results from the various new methods described in Section 3. 
Although dEJMH is based on the same methods as eEJMH, we observe that the runtime still grows with 
the size of the structuring element. This is because we work with the inverted input image X’^ H REC^ and 
the rectangle REG® depends on the size of the structuring element. Still eEJMH and dEJMH seem to 
provide a better alternative in most cases. 

8. Conclusions 

In this paper we developed new ideas to speed up the calculation of erosion and dilation. We also proposed 
fast algorithms to perform these operations with arbitrary structuring elements. Further we determined the 
runtime complexity of our algorithms which is given hj 0{\X\M) for erosion and by 0(|X'^nREC® \'M) for 
dilation, where X is the input image, B the structuring element, M denotes the number of runs of R, and 
REC^ is a rectangle of the size of X plus a border that is twice the size of B. Finally experiments confirmed 
that our algorithms provide huge speedup compared to some other well known implementations. 
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RLE Machado ^ OpenCV eEJMH 


dRLE dOpenCV dEJMH 



20 40 60 80 100 120 140 160 180 200 





Size of the structuring element (px) 


Size of the structuring element (px) 


Figure 2. Upper left: Erosion with a square-shaped structuring element. Upper right: Dilation with a square-shaped structuring 
element. Lower left: Erosion with a diamond-shaped structuring element. Lower right: Dilation with a diamond-shaped 
structuring element. 
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